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Abstract 

Motivated by the problem of channel estimation in wireless com- 
munications, we derive a reconstruction formula for pseudodifferential 
operators with a bandlimited symbol. This reconstruction formula 
uses the diagonal entries of the matrix of the pseudodifferential oper- 
ator with respect to a Gabor system. In addition, we prove several 
other uniqueness theorems that shed light on the relation between a 
pseudodifferential operator and its matrix with respect to a Gabor 
system. 

1 Introduction 

The mathematical formulation of orthogonal frequency division multiplexing 
(OFDM) in wireless communications uses several fundamental notions from 
time-frequency analysis. On the one hand, Gabor expansions are used to 
transform digital information into an analog signal. On the other hand, 
pseudodifferential operators are used to model the distortion of a signal by 
the physical channel. Inevitably the rigorous analysis of the communication 
system leads to new and interesting questions in time-frequency analysis 
that are quite relevant for communication engineering. 

In this paper we study a problem arising in channel estimation. Which infor- 
mation is required to determine the symbol of a pseudodifferential operator? 
How can an operator be reconstructed from such information? 



"This work was funded by the Austrian Science Fund (FWF) in project NFN SISE 
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To put the discussion on a firm basis, let us describe an extremely simplified 
model of signals and the transmission in wireless communications. See [161 
H3] for the mathematical models. Let n(z)f(t) = e 2mZ2t f{t — z\) denote 
the time-frequency shift by z = (z\, z<i) £ R 2 , and let A = aL x bZ denote 
a lattice in the time-frequency plane with lattice parameters a, b > 0. In 
orthogonal frequency division multiplexing (OFDM) a string of numbers c\, 
the "digital" information, is used as the coefficient sequence for a Gabor 
series of the form 

f = J2J2 Ckie 2mblt g(t - afc) = J] cavt(A) 5 , t € R . 
k£Z\l\<B AeA 

The pulse g is usually taken to be a characteristic function, but in nonsta- 
tionary environments pulses with better frequency concentration are prefer- 
able [2j[T71[18l[231[37] . The analog signal thus built is then transmitted from 
a sender to a receiver and distorted or transformed by physical processes. 
The second link between wireless communications and time-frequency anal- 
ysis is the description of the distortion of the signal / during the physical 
transmission. As a result of multipath propagation and of the Doppler effect, 
the received signal is a superposition of time-frequency shifts. Specifically, 
the received signal can be written as 

f(t) = J J £(77, u)tt(-u, r/)f{t) dudr] . 

Here a is the Fourier transform of a symbol a on IR 2 and is called the spread- 
ing function that indicates the amplitude of each occuring time-frequency 
shift. In the standard mathematical language the distortion / — > f is just 
the pseudodifferential operator (in the Kohn-Nirenberg calculus) with sym- 
bol cr, and is usually written as 

f(t) = a KN f(t) = J a(x,0e 2m ^f(0dC. 

For physical reasons the time delay and the Doppler shift must be bounded, 
and therefore the spreading function a has a compact support. Equivalently, 
the symbol a is bandlimited, i.e., analytic and of exponential type. From 
the perspective of analysis, such pseudodifferential operators are extremely 
special and are only the raw material for the study of difficult operators [21] . 
For wireless communications, pseudodifferential operators with bandlimited 
symbols are precisely the appropriate model. 
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At the receiver, the distorted analog signal / is analyzed by taking corre- 
lations with time- frequency shifts of the given pulse (or some other pulse). 
Thus the data to be analyzed are therefore the numbers 

yx = (f, Tr(X)g) = c fl (a KN 7T( f ,)g, tt(%) A <E A. (1) 
iueA 

The task of the engineer is now to recover and estimate the original data c\ 
from the received information y\. The central object here is the matrix H 
with entries 

H XlA =((a KN ^) g ,n(\)g}) x ^ A (2) 

of the pseudo differential operator with respect to the set of time-frequency 
shifts {7r(A)g : A G A}. In wireless communications this matrix is called 
the channel matrix. Its estimation and inversion are among the principal 
engineering tasks. 

A fundamental mathematical problem concerns the relation between the 
channel matrix and the symbol. This range of questions has been studied in 
time-frequency analysis, e.g., in [UlCE], yet the models and assumptions of 
wireless communications pose new and intriguing problems. An important 
objective is to recover or approximate the symbol a from partial information 
about the channel matrix; this is the problem of channel estimation. Usually, 
in real wireless communication systems, pilot tones are used to estimate some 
entries of the channel matrix on the diagonal [6,44J. The problem then is to 
recover the entire matrix ([2]) and subsequently to solve the system y = He. 
The engineering models lead to the mathematical question when and how 
the channel matrix is completely determined by its diagonal. For arbitrary 
operators, this question does not even make sense, but for operators with 
an analytic symbol, as we will see, one can recover the symbol completely 
from the diagonal of the channel matrix. Our first contribution is a precise 
reconstruction formula for the symbol from the diagonal of the channel ma- 
trix. In other words, the matrix is completely determined by its diagonal! 
The hypothesis that a is band-limited suggests a connection to the sampling 
theory of band-limited functions. Indeed, once this connection is established 
(which we do in Lemma [3 . 1 1 below) . one may apply results from the Shannon 
sampling theory and obtains the following reconstruction formula. 

Theorem 1.1. Let A = aZ d x VL d , g G S(R d ). If a G S'(R 2d ), suppa C 

Q £ = [— ^7 + £) 2a ~ £ ] rf x t — 15 + e ' lb ~ £ ] d f or some e > 0, then there exists 
a kernel K G S(R 2d ) such that 

a = ^2(a KN 7T(X)g,7r(X)g)T x K (3) 
aga 
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with convergence in S' . 



Theorem 11.11 provides a theoretical answer to a crucial point of channel 
estimation: how can the channel be estimated from (partial) knowledge of 
the diagonal of the channel matrix? 

We formulate several versions of this theorem that reflect various models 
and assumptions in wireless communications. In engineering it is usually 
assumed that the channel is a Hilbert-Schmidt operator and a G L 2 (M. d ). 
In this case, the series (|3j) converges in L 2 , and ([3]) is valid under weaker 
assumptions on the pulse g. In our opinion, however, the distributional 
version offers a better model for signal propagation since "point scatterers" 
correspond precisely to point measures in the symbol. Moreover, for the 
stable recovery of the coefficients c\ from ([I]), the channel matrix must be 
invertible, which is certainly not the case for a Hilbert-Schmidt operator. 
Theorem 1 1 . 1 1 can be interpreted as a result about operator identification and 
shares several aspects with the work of Pfander et al. [25^32] . In a series of 
papers, they investigate the question under what conditions a symbol a can 
be recovered from a single measurement o~ KN f on a suitable distribution /. 
Similarly to Theorem 11.11 their answer is expressed as a sampling theorem 
and is valid for bandlimited symbols. 

Our second contribution is the analysis of some common assumptions in the 
engineering community. A common assumption in the wireless communica- 
tions literature is that the channel matrix is diagonal so that its inversion 
becomes trivial. Equivalently, this means that the channel matrix is di- 
agonalized by the time-frequency shifts of a suitable function. There are 
numerous papers building on this assumption, e.g. [HIEKSZIEH] ■ We show 
that this assumption cannot withstand mathematical scrutiny. 
We prove that if the underlying Gabor system is a frame for L 2 (R d ), but 
not a basis, then the corresponding channel matrix H cannot be diagonal, 
unless the operator a KN is identically zero. We further prove that for a non- 
zero pseudodifferential operator with a bandlimited symbol and a Gaussian 
window, the channel matrix cannot vanish identically, quite independently 
of the spanning properties of the Gabor system. 

The paper is organized as follows: In Section [21 we summarize the mathe- 
matical preliminaries. In Section [3] we present the reconstruction formula for 
the symbol of the pseudodifferential operator, and in Section |3] we present 
the proposed uniqueness results. 
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2 Preliminaries 



We collect some concepts and definitions from time- frequency analysis. The 
precise details and proofs can be found in [11] or in |13| . 
The Fourier transform of a function / E L 2 (M. d ) is defined as 



ho = I f(x)e~ 2 ^- x dx, e e 



The two fundamental operators in time-frequency analysis are the transla- 
tion operators T x and the modulation operators defined by 

T x f(t) = f(t-x) and M^f(t) = e 2 ^ t f(t), t,x,^€R d . 

Their compositions are the time-frequency shift operators ir defined as 



n(z)f(t) = M,T x f(t) = e 2 ^-*/(i - x), for z = (x, £) G R 2d . 



A set of time- frequency shifts of a non-zero window function g £ _L 2 (R rf ) 
with respect to a lattice A = aZ rf x 6Z d C M 2d , a, b > 0, 

G{g,A) = {iT(X)g : A € A} 

is called a Gabor system. If there exist constants A, B > such that for all 

^ii/ii 2 <EKM A )2>i 2 < 5 imi 2 ' 

AeA 

then the set G(g,A) is called a Gabor frame with frame bounds A and i?. 
The short-time Fourier transform (STFT) of a function or distribution / 
with respect to a non-zero window g is defined as 



= / f(t)g{t-x)e- 2 ^ t dt 

JR d 

= (f,M 6 T x g) = {f,ir(z)g), 

for z = (x, £) G M M . 

The Rihaczek distribution of two functions f,g€ L 2 (M. d ) is defined as 



= /W(0e-^, (4) 

for x, £ 6 R . The Rihaczek distribution and the short-time Fourier trans- 
form are related in the following way: 

RUTg) = uv g f, (5) 
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where UF{£,x) = F{-x,£) and x,£ G R d . 

Let a be a (measureable) function or a tempered distribution on R 2d . The 
bilinear form 

(a KN f,g) = (a,R(g,f)), f,gGS(R d ), (6) 

defines a linear operator o~ KN from to S'(R d ). The operator a KN 

is a pseudodifferential operator in the Kohn-Nirenberg calculus with Kohn- 
Nirenberg symbol a. If a is a locally integrable function, then the Kohn- 
Nirenberg transform can also be written as 

a KN f(x)= [[ a(r ] ,u)M ri T^ u f(x)dud? ] . (7) 

In engineering this version is called the spreading representation of the pseu- 
dodifferential operator and a in (J7J) is known as the spreading function of 
o~ KN , since it describes how much the function / is "spread out" in time 
and frequency under the action of a KN . 



3 Reconstruction Formula 

We first deal with the question if and how the symbol of a pseudodifferential 
operator o~ KN can be reconstructed from the diagonal of the channel matrix. 
Let A = aZ d x bZ d be a lattice in R 2d , g G L 2 (R d ) a non-zero window 
function, and Q(g, A) be the corresponding Gabor system. The matrix H of 
a pseudodifferential operator a KN with respect to the Gabor system Q(g,A) 
is defined as follows: 

H Xll = (a KN 7r(fi)g,7r(X)g), X, € A. 

When a KN describes a wireless channel, the matrix H describes the action 
of a wireless channel on certain transmit pulses and is therefore called the 
channel matrix. 

We remark that for the solution of the linear equation (pQ) H must be in- 

vertible. In wireless communications it is costumary to assume that Q(g,A) 

is a Riesz basis for the generated subspace p4l[27ll28l l42] . 

We first derive an alternative expression for the diagonal entries of H in 

terms of the Rihaczek distribution of g. Let J 7 L 1 (R 2d ) denote the Fourier 

algebra on R 2d consisting of all functions on R 2d with integrable Fourier 

transform. 

We have the following well-known lemma. 
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Lemma 3.1. Let a G V 1 {^?' d \\ < p < oo, a compactly supported, g G 
nfX 1 ^). The diagonal entries of H can be written as follows 



H x ,x = (a KN ir(\)g,ir(\)g)=CT*R(g,gy(\), A G A, 



where f*(x) = f(—x). 

Proof. Under the given assumptions, a is infinitely differentiable and D a a is 
bounded for all multi-indices a. The standard theory of pseudo differential 
operators implies that a KN is bounded on L 2 (R d ) [1 lj.121] . Consequently 
the mapping A G M. 2d — > (a KN ir(X)g, Tr(X)g) is continuous, and the channel 
matrix is well-defined. 

Using the definition of the Rihaczek distribution ([4]), we get 

\\R(g,9)* Hi = <?lli = \\g\\i ■ \\g\\i < oo. 

Since a G L p (M. 2d ), the convolution a * R(g,g)* is well-defined in the V- 
sense. 

From the intertwining property of the Rihaczek distribution |15l Lemma 
4.2], we have 

R(7T(X)g,7r(X)g)(z)=R(g,g)(z-X). (8) 

Combining the definition of the Kohn-Nirenberg transform ([6]) and Equation 
([8]), we obtain 

(o- KN n(X)g,n(X)g) = (a, R(w(X)g, w(X)g)) 



a(z)R(ir(X)g,ir(X)g)(z)dz 



o-(z)R(g,g)(z - X)dz 

a(z)R(g,g)*(X- z)dz 
= a*R(g,g)*(X). (9) 
In general, ([9|) is valid for almost every A G M. 2d . Since 

supp T\ o- * R(g, g)*) = supp (a ■ R(g, g)*) C supp a , 



is compact, a*R(g,g)* is an analytic function, and therefore © is valid for 
every A G M? d . U 
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By combining the observation of Lemma l3.1l with the classical sampling the- 
orem for band-limited functions, we obtain a reconstruction formula for the 
symbol of a pseudodifferential operator from the diagonal of the channel ma- 
trix. In the formulation of the multivariate version of the Shannon-sampling 
theorem, we need the "sine" -function adapted to a lattice A = aL d x b7L d C 
R M , namely 

d . 2d . , 

TT smiraxj j-r smiroxj 

smc(i) = . 

- LJ - iraxi - LJ - irbxn 
j=i J j=d+l J 

Then every function a G L 2 (R 2d ) with supp / QQ = [-^, ±} d x 
possesses the cardinal series expansion 

/ = ^/(A)T A sinc 

AeA 

with convergence in L 2 and uniformly. For the general theory of Shannon 
sampling we refer to [5t[2TJ|[50] . 

The first reconstruction formula is stated for symbols a in L p (R 2d ), 1 < p < 
oo. 

Theorem 3.2. Let a G I^(R 2d ),l < p < oo, suppa C Q = \-^.,h] d x 



2a' 2a 

^, ±] d and g G L 1 ^) n J^X 1 ^). Choose ip G C™(R 2d ) such that tp = 1 



on Q, define K = T y uv g ) an ^ assume ^ ia ^ ^V g g does not vanish on 
suppip. Then the symbol a can be reconstructed from the diagonal entries 
H\ \ = (a KN Tr(\)g, 7r(A)<?) of the channel matrix via the modified cardinal 
series 

a = -— 1 ^2H X! xT x (sinc*K). (10) 

^ a > AeA 

The sum converges absolutely, uniformly and in V{R 2d ). If p = 1, the 
reconstruction formula is still valid, but the series converges only in L q (R 2d ) 
for q>\, but not in L 1 (R 2d ). 

Proof. We apply the multivariate version of the classical Shannon-Whittaker- 
Kotelnikov sampling theorem with the lattice A = a% d x VL d to the ban- 
dlimited function a * R(g, g)* . 

We recall from Lemma 13. II that H\ t \ = a * R(g,g)*(\) and write 

a * R(g, g )* = -L^J2 Hx,xT x smc, (11) 
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According to the V-theory of the cardinal series [HE] the sum converges 
absolutely, uniformly, and in ^(M 2 ^) for 1 < p < oo. This sampling ex- 
pansion holds pointwise also for p = 1, but the convergence is then only in 
L q (R 2d ) for q > 1. 

Since R(g,g)* G L 1 (M 2d ), formula © implies that UV g g = R[g~7g)* G 
T L 1 (M. 2d ) . Since UV g g / on supp</?, we apply the Wiener-Levy the- 
orem, see |34} Theorem 1.3.1], to conclude that there exists a function 
ip G J r L 1 (R 2d ) such that ip = on supp p. Since <p G C™(R 2d ), cp is also 

UV g g 

in TL x (R? d ) and we have ipi\) = ^f=- G J r X 1 (M 2d ). Thus we conclude that 
Now we claim that 

a*R(g,g)* *K = a. (12) 
To prove this, we take the Fourier transform of (112D and obtain 



& ■ R(g, g)* ■ K = a ■ UV g g ■ UV g g 1 • ip = a, 

because ip = 1 on supp a. 

Finally we combine (fTT]) and (fl~2j) to compute 

a = a * R{g, gf * K = ^ H x<x T x ( S mc * K). (13) 

Since we convolve a * R(g,g)* G -^(M 2 ^) with K G i 1 (M 2d ), the series in 
(|13p inherits the convergence properties from (lll|) . I 

Similar reconstruction formulas are valid when one or several side-diagonals 
of the channel matrix H are known [31} Chpt.4]. In the engineering practice 
the diagonal entries are estimated only on a sublattice A p C A by means of 
pilot symbols [6]l44j . In this case one applies Theorem 13.21 to a sublattice of 
the full lattice. 

In the case p = 2 we can weaken the assumptions on g considerably to g G 
L 2 (M. d ). This case is important because it treats the (unjustified) assumption 
of the engineering community that wireless channels are Hilbert-Schmidt 
operators [251(33] . In addition, it covers the rectangular window g = X[ a ,/3] d 
corresponding to OFDM without pulse-shaping [6j ll0t l4"4"]. 

Proposition 3.3. With the notation of Theorem \3.2\ assume that a G 
L 2 (R 2d ) and g G L 2 (R d ). Then a can be reconstructed from (Hx t \)\eA 
by 

a = -j-Y J Hx,\T x (sinc*K). (14) 
^ a ' AeA 
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with convergence in L 2 (M. 2d ) and uniform convergence. 



Proof. The proof requires only a minor modification. Since R(g,g)* = lAV g g 

is bounded, the product a ■ R(g,g)* is in L 2 (M. 2d ) with support in Q. Thus 
a * R(g,g)* G L 2 (M. 2d ) is bandlimited and the sampling reconstruction (jlip 
holds with uniform convergence and convergence in L 2 (R 2d ). 
Finally, since UV g g does not vanish on supp <p by assumption, the multiplier 
K = tp ■ UV g g is bounded, and therefore the operator F h- > F * K is 
bounded on L 2 (R 2d ). Consequently, the deconvolution formulas (fT2j) and 
(I13p are well-defined on L 2 (IR 2d ) and the reconstruction (I14D follows. I 



Finally we formulate a distributional version of the reconstruction theorem. 
This version is not just for the sake of mathematical generalization, but is 
necessary for the accurate modelling of physical channels. For example, a 
single point scatterer with time delay r and Doppler shift v has the point 
measure S( T>V ) as its spreading function. A typical spreading function is 
usually written as a distributional part plus a random component [26]. By 
adapting the hypothesis of Theorem 13.21 we obtain the following statement. 

Proposition 3.4. Let a G S'(R 2d ), suppa C Q £ = [-i + e, ± - e] d x 

[~3B + e > ^5 - £ ] d /° r some e > and 9 G «5(M d ). CTioose <p G ^(IR 2 ^) suc/i 
i/iai = 1 on Q £ , suppip C [— x [— and assume that UV g g 
does not vanish on supp (p. Then the symbol a can be reconstructed from the 
diagonal entries H\\ = (o~ KN ir(\)g,ir(\)g) of the channel matrix via the 
modified cardinal series 

a = J—y Hx xTx T- 1 ( , (15) 
^ V UV g g J 

with distributional convergence. 

Proof. If g G S(R d ), then R(g,g) G S(R 2d ) and thus a * R(g,g)* G S'(R 2d ) 
with supp a * R (g, g)* C Q e . The distributional version of the sampling 
theorem [1] now yields that 

a * R(g, 9 y = ^(a * R(g, g)*)(\)T\J r ~ 1 (p. 

with distributional convergence. Since cr*R(g,g)* is an entire function of at 
most polynomial growth on M? d (by the theorem of Paley- Wiener [35]) the 
pointwise evaluations are well-defined. Likewise, since a KN is continuous 
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from 5(R d ) to the mapping A i-> (a KN Tr(X)g,Tr(X)g) is continuous, 

therefore, as in Lemma [3TTI H\ t \ = (a * R(g, g)*)(X) for A £ A. 

To conclude, we observe that ip • UV g g 1 is in 5(IR 2a! ) and thus also K = 

(^p ■ UV g g ^ G 5(M 2d ). Consequently, the deconvolution formulas (fT2|) 

and (|13p make sense in 5'(IR M ), and the reconstruction formula is proved. 
■ 

Proposition 13.41 will also be relevant for the numerical implementation of 
the reconstruction formula. By assuming a slightly smaller spectrum Q e , 
the expanding kernel J 7-1 (^p -UV g g ^ is in 5(M 2rf ) and decays rapidly. If 

a G LP{E? d ) instead of 5'(R d ), then the expansion is localized and converges 
rapidly. 



4 Uniqueness Results 

The reconstruction results of the previous section imply that a pseudodif- 
ferential operator with a bandlimited symbol is uniquely determined by the 
diagonal of the channel matrix. In this section, we prove further uniqueness 
results that illustrate the relation between a pseudodifferential operator and 
the corresponding channel matrix under various assumptions on the Gabor 
system and the symbol. For notational simplicity, we now work in dimension 
d=l. 

Let A = aL x VL be a lattice, let g G X 2 (M) and Q(g,A) the corresponding 
Gabor system in _L 2 (R). We assume that 

H Xll = (a KN 7r(^g,7r(X)g) =0, VA, /i€ A, (16) 

Of course, if G{g-, A) spans L 2 (M. d ), then obviously a KN = 0. Under the basic 
assumptions of wireless communications (bandlimited symbol and Gabor 
Riesz sequence) the conclusion is not so obvious. 

Before we state the next theorem, we rewrite the general entries of the chan- 
nel matrix. We write A = (Ai, A2), fJ. = ([J,i, ^2) 6 From the definition of 
a , see ([6]), we have 

(a KN n^)g,7r(X)g) = (a,UV^ )g ir(X)g). 

Using the covariance property of the STFT (e.g., [131 Ch. 3]), we compute 

V nWg n(X)g(-x,0 = e- 2 ^- x ^e 2 ^- x - x ^T x ^V g g(-x,0 
= e^ iX ^e- 2m ^M { ^_ Xl) T^^V g g{-x^) . 
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Writing G = UV g g, we conclude that (fTB]) is equivalent to 

(a,M (M2i _ Al) T A _ M G)=0, VA, /i G A, (17) 

Here // 2 G frZ and Ai G aZ, so the modulations are taken with respect to the 
lattice A' = VL x aZ. Combining (|16p and (|17p . we obtain 

(<7, M^T x G) = 0, VA G A, n G A' = 6Z x aZ. (18) 

The first uniqueness theorem treats the case of Gaussian Gabor systems. 

Theorem 4.1. Let a G 5'(M 2 ), ip(x) = e -7 ™ 2 be the Gaussian function and 
A = aZ x 6Z an arbitrary lattice. If a is compactly supported and 

H\[i = (a KN Tr(ti)<p, TT(X)ip) = 0, VA, ii G A, (19) 

then a KN is identically zero. 

Proof. The STFT of the Gaussian function is given by |13t Lemma 1.5.2] 

G(£,x)=UV v <p(S,x) = V>(-x,£) = -L e -|« 2 -f- 2 +^. 

Setting n = Qui, /U 2 ) G A', A = (Ai, A2) G A, z = (£, x) G M 2 , we obtain that 
M^T X G = e 2 ™ z ^G{z - A) 

V2 



The exponent in (I20p can be written as 

27ri(^i + xp 2 ) + 7ri(£ - Ai)(s - A 2 ) - |(C - Ai) 2 - |(x - X 2 pl) 

7T 2 7r 2 ^2 ^2 

= 7riAiA 2 - -Ax - -A 2 - -£ - -x +m£x- 
— 7rz£A 2 — TrixAl + vr^Ai + 7rxA 2 + 2m{^pLi + xp 2 ). 

Next we define 

^a(£ ^) = e" 7 ™^ 2- 7 ™ :E ^i+ 7r £^i+ 7r:c ^2 A G A 

and 

p = a ■ G. (22) 

Since G is in 5(M 2 ), p is a tempered distribution with compact support 
in a set S C ]R 2 . Furthermore, since Fa is infinitely differentiable on any 
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bounded open set, F x p is again in 5'(M 2 ) with compact support in S. With 
this notation and k\ = e 7riAlA2_ 2 A i~i A 2 ^ 0, the assumption (fT9j) can be 
recast as 

(<7,M M T A G) = k x (a,M lx GF x ) = k x (p,M f ,F x ) 

= k x (p-F x )\p) = (23) 

for all A E A and for all \x € A'. 

In the remainder of the proof we will apply the Poisson summation formula 
for compactly supported distributions [121 Corollary 8.5.1] for each A G A. 
By restricting the periodization to a specific open set, we will derive a system 
of equations for the restrictions of p. Then we will use the invertibility of a 
Vandermonde matrix for some selected A's to conclude that all restrictions 
of p vanish. 

By the Poisson summation formula for compactly supported distributions, 
(|23j) is equivalent to 



&\ := E T »(P F *) = °> VAeA ' 



where A° = yL x is the adjoint lattice of A and thus the dual lattice of 
A'. Since S is compact, there is a positive integer L such that 

f-L + 1 L-l\ f-L + l L-l 

b = supp p C K 



We define the open rectangle Q = (— \, \) x (-±, ±) and Q jik = (|, | ) + Q, 
for — L < j,k < L. Clearly Q contains a period of £P X and 

-i+l<i,fe<i-l 

Next we consider the restriction of the periodization S? x to Q, 

&x\Q= E ^(pF A )| Q = 0. (25) 

Since i? contains supp p, (f25j) reduces to 

L L L L 

^A|Q= 2 E r fi 4)^)10= E E T (^)PlQ- T (i^) i? A|Q = 
j=—L k=—L j=—L k=—L 

(26) 
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for all A 6 A. In (|26p the sum is from — L to L, because this covers all shifts 
of Q which intersect supp p according to (f24"|) . These identities (as well as 
the subsequent derivations) are to be understood in the weak sense. For 
instance, (|26p means that for all ip G 5(M 2 ) with supp<£> C Q 

L L 

(&x\ Q ,<p)= E ^ T (|,!) /9 " T (i,!) i?A '^ =0 - 

j=-Lk=-L " 

Next we define 

Pj,k ■= T ( i k) p\Q . 

v b 1 a ' 

Since 

F x (Z-j/b,x- k/a) = F x ^,x)F x (-j/b,-k/a) (27) 

and F\(^,x) ^ for all we may divide (|26j) by the exponentials in 
x and £ and obtain 

L L 

E E ft'.* F x(-j/b,~k/a) 

j=-Lk=-L 

= E E ^-fe e^^-^'D^M-^!) = 

for all A = (aZi , 6Z2) £ A. By splitting the sum over k into an even and an 
odd part, we obtain 

E E pi* (e-^y>(^i) + 

j=~L k even 
j=—L k odd 

This identity holds for all h, l 2 G Z. We note that the (4L + 2) x (4L + 2)- 
matrix 1/ with entries Vi ir = z 1 ^ , r = 1, . . . , 4L + 2, li = 0, 1, . . . , 4L + 1, 
where z r = e ^-L-i)a/b f or r = i ; . . . , 2 L + 1 and z r = - e ™{r-ZL-2)a/b 
for r = 2L + 2, . . . , 4L + 2 is a Vandermonde matrix based on the 4L + 2 
distinct nodes ±e 7r - ?a ' /b . Likewise the (2L + 1) x (2L + l)-matrix W with 
entries W; 2 fc = e - l ^ kb / a } k = —L, . . . , L, l 2 = 0, . . . , 2L is an invertible Van- 
dermonde matrix. 

We conclude that, for every l 2 and for every j = —L,...,L, 

E Pi,k e 12 ^' - ^-) = (28) 

fc even 
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and 

E Pj,k e 12 ^'-**-) = 0. (29) 

k odd 

We divide the equations (f28|) and (f29|) by e /27nj and add them to obtain 

E Pi.* (e"^) 2 =0 (30) 

fc=-L 

for every I2 and for every j = —L,...,L. Since the coefficient matrix W is 
invertible, we conclude that pj & = for — L < j,k < L. 
Using [201 Theorem 2.2.1], we conclude that p = 0. Since p = oG by (|22|) 
and G is a Gaussian, we arrive at <r = 0, as was to be shown. 



Idea of an alternative proof of Theorem 14.11 

We sketch an alternative proof where the argument is based on the fact 
that the Gaussian g is a (strictly) totally positive function. This means 
that for two arbitrary sequences of real numbers x\ < X2 < ■ ■ ■ < x n and 
Vl < U2 < ■ ■ ■ < Vn the matrix 

{g{xj - Vk),j,k = 1, . . . ,n) 

has a strictly positive determinant and is thus invertible. Sec [39 , 40J for the 
fundamental properties of totally positive functions. 

It is profitable to use the Weyl calculus of pseudodifferential operators, which 
is formulated by means of the (cross-) Wigner distribution of /, g G £ 2 (M) 

W(f,g)(x,0= [ f(x+±)g(x-l)e~ 2m Z 4 dt. 

JR 

The Wigner distribution satisfies the following covariance property ( [11] 
or [El Prop. 4.3.2c]): 

W(n(p)n(X)fM^)9) = M A T A W '(/ ', g) , 

where p = (p2, —pi) for p = (pi,p2) £ The Wigner distribution of the 
Gaussian cp(t) = e~ wt is the Gaussian 

W(cp,<p)(x,£) = V2e^ x2+e ^ 2 

with Fourier transform <!?(£, x) = W(<p,<p)(£,x) = 2~i e - 27r (? 2 + a ' 2 ). Now 
assume that a € X 1 (IR 2 ) and set f(£,x) = e nlx ^a(^,x), then we have 

{a KN ir(\)g, n(p)n(X)g) = (r, ^(7r( M )7r(A) ff , 7r(A)<?> (31) 
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(See [IT] or [1 3|. Chpt. 14.3] for the transition between the Kohn-Nirenberg 
calculus and the Weyl calculus.) Consequently, (a KN Tr(i')(p,7r(X)ip) = for 
all A, v G A holds, if and only if 

(r, W(7r(MMA)p, tt(A) V ) = e™ 2 (r, M A T A tf%, p)) 

= e ^ lM2e -2^A-A ^ . t a $)^(A) = (32) 

for all A,/i G A, or equivalently for all A G A and /i € A'. By the Poisson 
summation formula f|32|) is equivalent to 

^eA° j,k=—L 

for almost all x,£ G [0, 1/a] x [0, 1/6] and for all /i = (bfa, —al\) G A'. Since 
Gaussian functions are totally positive, the matrix with entries e~ 7T ^ + i~ bl ' 2 ' ) I 
\j\ < L, h = 0, . . . , 2L is invertible, and likewise the matrix e - 7r ( !E +i'~ oZl ) )/ 2 , 
< ^2 = 0, . . . 2L is invertible. From this we conclude that t(£ + ^, x + 
I ) = for all j, k and almost all x, £. Consequently f = almost everywhere 
and thus a = 0. ■ 



Finally, we investigate an important assumption made in wireless commu- 
nications. To facilitate the inversion of the channel matrix in ([I]), it is com- 
monly assumed that the channel matrix is a diagonal matrix 8. !). 271138] . 
Clearly, if the Gabor system G(g, A) is a (Riesz) basis for L 2 (M. d ), then there 
is a bijection between operators and channel matrices, and thus there exist 
diagonal channel matrices. 

However, if the Gabor system is a frame, then this assumption is never 
satisfied. 

Theorem 4.2. Let a KN be a bounded operator on L 2 (R d ), A = aZ d x VL d 
be a lattice with ab < 1 and g G L 2 (M. d ) such that Q(g,A) is a frame for 
L 2 {R d ). If 

(a KN 7r(n)gM^9) =d„d Xlx , VA, /i G A, (33) 
then a KN is identically zero. 

Proof. It follows from f|33f) that a KN Tr(fj,)g is orthogonal to the linear span 
of the set {ir(\)g : A G A,A / /i}, V/x G A. That is 

a KN n(fj,)g _L span{7r(A)# : A G A, A / fj,}, V/x G A. 
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Since {7r(A)<7 : A £ A} is a frame for L 2 (R d ), according to Lemma IX], 
there are two possibilities for the set {7r(X)g : A G A, A / fi}. Either 

• {7r(A)(? : A G A, A / //} is incomplete in L 2 (R d ), or 

• {7r(A)(? : A G A, A /i} is again a frame for L 2 (M d ). 

If the set {-7r(A)g : A £ A,A / /i} is incomplete in L 2 (M. d ), then the set 
{7r(A)(? : A G A} has to be a Riesz basis for L 2 (M. d ). Otherwise it would not 
be possible to have an incomplete set after removing one single element. It 
follows from the density and duality theory of frames [5, 13,22] that in this 
case ab = 1. This contradicts the assumption ab < 1. 

Therefore the set {ir(\)g : A G A, A / fi} is again a frame for L 2 (R d ). This 
implies that the orthogonal complement of {ir(X)g : A G A, A / /i} has to be 
equal to the set {0}. That is, 

a KN Tv(fx)g _L span{7r(A)# : A G A, A ^ fj,}, V/i G A, 

is only possible if a KNr K{n)g = 0, \/fi G A and therefore a KN = 0. ■ 

It seems that some fundamental algorithms of wireless communications in 
nonstationary environments are based on the incorrect assumption that the 
channel matrix is diagonal. Nevertheless the intuition of the communication 
engineers is perfectly correct and can be supported by rigorous mathematical 
results. Indeed, if the symbol a is smooth and the pulse g of the Gabor 
system possesses a minimal amount of time-frequency concentration, then 
the channel matrix decays rapidly off its diagonal |141ll6j. Hence, from a 
numerical point of view, the channel matrix can be approximated well by 
a diagonal matrix. This idea was used for improved equalization methods 

in [rang. 
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